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In the present work, flat FRW model of the universe is considered to be an isolated 
open thermodynamical system where non-equilibrium prescription has been studied using 
the mechanism of particle creation. In the perspective of recent observational evidences, 
the matter distribution in the universe is assumed to be dominated by dark matter and 
dark energy. The dark matter is chosen as dust while for dark energy, the following choices 
are considered: (i) Perfect fluid with constant equation of state and (ii) Holographic dark 
energy. In both the cases, the validity of generalized second law of thermodynamics (GSLT) 
which states that the total entropy of the fluid as well as that of the horizon should not 
decrease with the evolution of the universe, has been examined graphically for universe 
bounded by the event horizon. It is found that GSLT holds in both the cases with some 
restrictions on the interacting coupling parameter. 
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I. INTRODUCTION 


Cosmological models using irreversible thermodynamics along with Einstein’s general relativity 
(GR) have recently gained great interest among cosmologists. With this new perspective, it is possible 
to analyze different cosmological models which have blockages in equilibrium thermodynamic set-up. 
Eling, Guedens and Jacobson [l] first introduced irreversible thermodynamics in cosmology after he 
failed to reproduce Einstein’s field equations from the first law of thermodynamics in f(R )-gravity and 
concluded that a non-equilibrium thermodynamic treatment is essential for curvature correction to 
entropy. They added an extra entropy term diS called entropy production term in the entropy balance 
equation dS = ^ + diS, where they have argued diS as the bulk viscosity production term which can 
be determined by imposing energy-momentum conservation. In general, the entropy balance relation 
in non equilibrium thermodynamics is of the form dS = d e S + diS, where d e S is the rate of entropy 
exchange with the surroundings while diS (> 0) comes from the process occuring inside the system. 
In particular, diS is zero for a reversible process and positive for an irreversible process. In cosmology, 
diS has no clear interpretation as it depends on the internal production process. This term has been 
identified as a dissipative effect by exploiting a fluid dynamics description of the local causal horizon 
kinematics while considering irreversible thermodynamics. Subsequently, a lot of work have been done 
in order to interpret the entropy production term in various ways. In 2009, G. Chirco and S. Liberati 
[2], have shown that the dissipative character leading to non-equilibrium spacetime thermodynamics 
is actually related (both in GR as well as in /(i?)—gravity) to non-local heat fluxes associated with 
the purely gravitational/internal degrees of freedom of the theory. In particular, in the case of GR 
they showed that the internal entropy production term is identical to the so called tidal heating term 


of Hawking-Hartle |3(. Similarly, for the case of /(R)-gravity, the dissipative effects can be associated 
with the generalization of this term plus a scalar contribution whose presence is clearly justified within 
the scalar-tensor representation of the theory. 


Two most profound mysterious components of matter present in the universe are dark energy (DE) 
and dark matter (DM). According to recent observational data obtained from Type la Supernova, 
our Universe is now experiencing an accelerated expansion. This late time cosmic acceleration is 
interpreted by the existence of DE which is believed to have a huge negative pressure and therefore 
it is gravitationally repulsive. A common and natural candidate for DE is the cosmological constant 
but it is not widely used in the literature due to its limitations and drawbacks. Several alternative 
dynamical DE models have been proposed but the nature of DE (except that it has a huge negative 
pressure) is still a mystery. On the other hand, the behaviour of galatic rotation curves and the mass 
discrepancy in the cluster of galaxies are interpreted by the existence of DM which unlike DE, is 
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gravitationally attractive in nature. As we know that DE is homogeneously distributed while DM is 
inhomogenously distributed, occuring in clumps around ordinary matter. Interaction between these 
two dark components is expected to be weak or even negligible. Still, the effect of interaction between 
them on Universal dynamics cannot be ruled out. In fact, recently it has been argued by Wang et al. 4] 
that an appropriate interaction between DE and DM plays a significant role in perturbation dynamics 
and affect the lowest multipole of the CMB sprectum. Harko and Lobo {5] have also investigated the 
interaction between DE and DM in the context of irreversible thermodynamics of open systems with 
matter creation/annihilation. In this connection, it should be mentioned that the thermodynamics of 
irreversible processes was introduced in cosmology by Prigogine et al. Jf|. 

In our present work, we have investigated non-equilibrium thermodynamics of the Universe bounded 
by event horizon. We have examined the entropy production term by considering the non-equilibrium 
process of particle (DE/DM) creation/annihilation as the internal production process. Here we have 
considered an isentropic system consisting of two interacting dark fluids. Karami and Gaffari [7] 
have examined the validity of the generalized second law of thermodynamics (GSLT) in a non-flat 
FRW universe containing interacting DE with cold DM in irreversible thermodynamical context and 
they have shown that GSLT is satisfied for certain range of energy-transfer constants. Here we have 
employed the theory of particle creation for describing interaction between DE and DM and examined 
the validity of GSLT for Universe bounded by the event horizon. 


II. PARTICLE CREATION MECHANISM IN THE TWO DARK SPECIES: A UNIFIED 

MODEL 


Suppose a closed thermodynamical system with N particles have internal energy E. The conser¬ 
vation of this internal energy gives rise to the first law of thermodynamics as 


dE = dQ — pdV , 


( 1 ) 


where p and V are the usual thermodynamic pressure and comoving volume and dQ stands for the 
amount of heat received by the system in time dt. Equivalently, we have the Gibbs equation as 


Tds = dq = d ( — ^ + pd (— 
\nj \n 


(2) 


where ’s’ is the entropy per particle, p = y is the energy density, n = y is the particle number density 
and dq = ^ is the heat per unit particle. It should be noted that the above Gibbs equation is true for 
an open thermodynamical system for which the number of fluid particles is not conserved 7^ 0 
This can be expressed mathematically as 


h + On = nT. 


( 3 ) 
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Here N 11 = nu^ is the particle flow vector, t/ 1 is the particle four-velocity, 0 = id), is the fluid expansion, 
r stands for the rate of change of the particle number in a comoving volume V and n = n^u^ by 
notation. T effectively behaves as a bulk viscous pressure which causes the thermodynamics to be 
non-equilibrium. T > 0 indicates creation of particles while annihilation of particles correspond to 

r < o. 

In the present work, we consider an open thermodynamical system to be flat FRW model of the 
Universe with particle creation characterizing the non-equilibrium phenomenon. Suppose that the 
Universe consists of two dark fluids namely DM and DE. The Einstein field equations are (choosing 
87tG = 1 = c) 

3-ff = Pt = Pm T Pd (4) 

and 

2H + 3 H 2 = —(jpm + n m ) — ( pd + U d ) (5) 

and the energy conservation equations are given by 

Pm T 3i7(p m + Pm) = 3-HTl m (6) 

and 

p d + 3H(p d + p d ) = -3HU d , (7) 

where (p m ,p m ,H m ) are respectively energy density, thermodynamic pressure and dissipative (bulk) 
pressure of DM and (p d ,p d ,H d ) are those of DE. The particle number conservation equations are 
modified as 


h-m T 3 Hn m — r m'li'T 


( 8 ) 


and 


h d + 3 Hn d = —T d n d , 


(9) 


where ( n m ,T m ) are the number density and particle creation rate for DM while ( n d ,T d ) are those for 


DE. We have assumed that T m > 0 and T d < 0 so that DM partic’ 


es are created while DE particles 


are destroyed, as predicted by the second law of thermodynamics js ]. If for simplicity, we assume the 
thermodynamical system to be isentropic (i.e., entropy per particle is conserved), then the dissipative 
pressures and the particle creation rates are related as [9] 


n m — q {Pm T Pm) 


(10) 
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and 

n d = ^j{pd + Pd)- (11) 

However, if the combined two-fluid is considered as a single fluid with energy density pt, thermody¬ 
namic pressure pt = p m + Pd and dissipative pressure = H m + n^, then combining (6) and (7), we 
have the usual conservation equation 

Pt + SH(pt + Pt) = ~ 3HHt. ( 12 ) 

Also, combining (8) and (9) we can write 

n t + 3Hn t = T t n t , (13) 


where nt = n m + rid is the total particle number density and T* is the particle creation rate of the 
single fluid. Further, if we assume the resulting single fluid to be isentropic, then we have 


At = —g(Pt + Pa¬ 


using (10) and (11) in (14), we obtain 


r* = 


r m{Pm 4" Pm) ^ d(Pd T Pd) 


(14) 


(15) 


(pt + Pt) 

The sign of T t indicates whether there is a creation or an annihilation of particles of the resulting 
single fluid. It should be noted that if H m = —< 0, then equations (6) and (7) imply interacting 
dark fluids with interaction term Q = 3HUd and the single fluid corresponds to a closed system as 

r t = o = n t . 


III. THERMODYNAMICAL ANALYSIS OF ENERGY TRANSFER BETWEEN THE 
DARK SECTORS OF THE UNIVERSE DUE TO PARTICLE CREATION 


Due to the observational evidences of the present late-time acceleration, the universe is currently 
dominated by a dark fluid system (DM+DE). The energy conservation relations for the subsystems are 
given by equations (6) and (7) respectively. Due to the particle creation mechanism, there is an energy 
transfer between the dark species and as a result, the two subsystems have different temperatures 
and thermodynamics of irreversible process comes into the picture. Thus, starting from the Euler’s 
thermodynamical relation nTs = p + p and using the energy conservation equations (6) and (7) as 
well as the modified particle number conservation relations (8) and (9) for the two dark subsystems, 
the evolution equations for the temperatures of the individual dark fluid elements are given by 


T, 


u=- 3 n“""+ifj+ 1+ .. 


OJn 


(16) 



6 


and 


where 


and 




\ dj d 

T) i + u d ' 


ef f m / -\ . 

q ( 1 H” Wn 


(17) 


(18) 


— u d + -^(1 + Wd). 


(19) 


Here, w m (0 < uj m < 1) and (< — |) are the equations of state parameters for DM and DE 
respectively and T m , T d are respectively the temperature of the two dark components. As for simplicity, 
the thermodynamical process is assumed to be adiabatic, so the effective bulk pressures H m and 
are respectively related to the corresponding particle creation rates by Eqs. (10) and (11). Now 
integrating Eqs. (16) and (17), we have the temperatures of the two dark sectors as 


T m = T 0 (l + u m )exp 


-3 


£ ( </f+ it) t 


( 20 ) 


and 


T d = T 0 (l + oj d )exp 


-3 


'a o 


u e J f - 


da 


( 21 ) 


where To is the common temperature of the two subsystems in equilibrium configuration and ao is the 
value of the scale factor in the equilibrium state. It should be noted that at the very early stages of 
the evolution of the universe, T m > T d and with the evolution of the universe, both the subsystems 
approach the equilibrium configuration (a = ao) with common temperature Tq. Then for a > ao,the 
thermodynamical equilibrium is violated due to the continuous flow of DE to DM and as a result we 
have T m < T$ < T d . In the present case, we are considering universe bounded by the event horizon 
as an isolated thermodynamical system and so in 
considered as the (modified) Hawking temperature 
where Re is the radius of the event horizon for the FRW model. Further, for thermodynamical 
analysis, we denote ’S TO ’ and , S d as the entropies of the two dark sectors of the present isolated 
system while Se stands for the entropy of the bounding event horizon. Then, from the first law (i.e., 
Clausius relation) for individual subsystems, we have 


the equilibrium state the temperature To may be 
10 , 111 on the event horizon, i.e., Tq = H 0 Re | a=a rM 


dS n 


dQr 


dE n 

dt 


T Pm' 


dV 

dt 


dt 


dt 


( 22 ) 
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and 


dSd = dQd = dEd dV 
dt dt dt ^ d dt 


(23) 


In the above expressions, E m = p m V and Ed = pdV are the energies of the two subsystems while Q m 
and Qd are the amount of heat of the two dark sectors and V = is the volume of the universe 

bounded by the two subsystems. Due to the isolated nature of the whole system, the heat flow across 
the horizon (Qh) is balanced by the heat flow through the two dark components, i.e., 


Qh 


(Qm Qd'j ■ 


(24) 


However, in view of non-equilibrium thermodynamical prescription one has to take care of the contri¬ 
butions from irreversible fluxes of energy transfer in the expression for the total entropy variation as 


IQ 


dSr 

dt 


dS m dSd_ dS E 
dt dt dt 


AdQdQd AhQhQhi 


(25) 


where Ad and Ay, are the energy transfer constants between DE and DM within the universe and 
between the universe and the horizon respectively. 

Here our universe is described by a 4—dimensional flat FRW model bounded by an event horizon. 
Normally in FRW spacetime, event horizon does not exist. It is relevant only when the universe is 
in an accelerating phase. So in the perspective of the present accelerating phase of the universe, it is 
reasonable to consider universe bounded by the event horizon as an open thermodynamical system. 
For event horizon, the area radius has the expression 

r°° dt 

R E = a - (26) 

Jt a 

(Note that the above improper integral exists when the universe is in an accelerating phase) and 
variation oi R E with respect to cosmic time t gives 


R e = HR e — 1. 


(27) 


Also, let us consider Bekenstein entropy-area relation on the surface area of the event horizon, i.e., 


S E = ttR 2 e . 


(28) 


Now we shall give two examples of DE model and discuss the above process in respective scenarios. 










A. Interacting dark energy with constant equation of state 


Here we shall consider DE as a perfect fluid with an equation of state p d = oj d p d where oj d is a 
constant such that uj d < —\ and DM to be in the form of pressureless dust, i.e., p m = 0. We choose 
dissipative pressures of DM and DE as 


n m — n d — a pd 


(29) 


so that the particle creation rates for DM and DE become (from Eq. (10)) 

r m = 3HX— , T d = 

Pm 1 + 


respectively. The temperatures of these two dark sectors read as 


T m = T 0 and T d = T 0 (l + u d ) (—) 

\Q> 0 / 


- 3 ( A +i^+^) 


Then, one gets (assuming 87 t = 1 = G) 


(30) 


(31) 


Q m = -(HR E )\\n d HR E + Q d - 1 ), 


Qd = --( HR E ) 2 n d (XHR E + u d + 1 ), 


(32) 

(33) 


and 


Q h = l(HR E ) 2 (l + u d Sl d ). (34) 

Differentiating (33) and (34), we have 

Qd = —H^R E Q d [XHR e (1 + qHR E ) + (A HR E + u> d + 1) {2(1 + qHR E ) + 3 HR E (uj d + a) d Q d + A + 2)}] 

(35) 

and 


Qh — —~H~R e [3HR E (jj d Qd(ud + to d Q d + A + 2) + 2(1 + u) d ^d)( 1 + qHR E )\ (36) 

respectively, where the deceleration parameter q can be evaluated as 

q = ~ 1 = ^(1 + 3uj d^d)- (37) 

Also, 

S E = i R e (HR e - 1). (38) 
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FIG. 1. The plots show the evolution of -Ar- against A for different values of ui d in Quintesssence (oj d = 
—0.4, —0.8), Phantom {u d = —1.2, —1.4) and Normal fluid {u d = §) eras. 


Thus, from Eq. (25), we get 


rlQE o 

T = —(HR E ) 2 (\n d HR E + n d -1) - 


dt 


2T n 

1 


27o (1 + aid) \aoJ 


\ 3(A+a; d + ( 1+ A „ d )) 


(HR E ) 2 n d {XHR E + u d + l) 
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+ -R e (HR e - 1) + -A d H 4 R E Q d (A HR e + u d + 1) XH 2 R E {1 + qHR E ) + (A HR E + u d + l) 


x {2(1 + qHR E ) + 3 HR E (uj d + uj d £l d + X + 2)} + ~AhH^R? E ( 1 + u> d Q d ) 
x (^d + u d £l d + A + 2) + 2(1 + u d Q d )(l + qHR E ) 


3 HR E uj d Q d 


(39) 


Due to a complicated expression of its evolution against A has been plotted in Fig. 1 for different 
choices of w d in quintessence, phantom and normal fluid eras. From the figure, we see that if the DE 
is either not close to ACDM or in phantom domain, then GSLT holds for all values of A while for DE 
in phantom domain, GSLT does not hold for very small values of A. Also, for w d = —1.4, the total 
entropy variation shows some peculiar fluctuations. 


B. Interacting Holographic Dark Energy Model 


The holographic dark energy density satisfies the relation p d = '-jj 13|, 14], where c is a constant 


5,Q 


to be determined from observations and L is an IR cut-off in units of 87 t = 1 = G. Li 


14] has shown 


that if we choose L as the radius of the future event horizon (R E ), only then we can obtain the correct 
equation of state of DE and the present accelerating universe. Moreover, interaction models have been 
shown to support [l5, 1^| the observationally measured phantom equation of state and are also favored 
by observational data from the cosmic microwave background (CMB) [l?]] and the matter distribution 
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at large scales [lg|. Thus, with the choice L = Re, we can write 


Re = 


(40) 


Hy/n? 

where VL d = Xjp is the DE density parameter. As in case A, here also we choose the dissipative 
pressures as — II m = 11^ = Xp = A p d , A being a small positive dimensionless parameter. Consequently, 
the particle creation rates for DM and DE take the forms T m = 3HX^j and T^ = — respectively. 


The density parameter evolves as [4j, 


m 


nl, = n n 




(i - n d ) - 3\n d 


(41) 


while the variable equation of state for holographic DE is given by 

1 2 /— 

Wd = — - — — — A. 

3 3c 


a 


(42) 


Note that in Eq. (41), the prime refers to derivative with respect to x = In a. Here the temperatures 
of DM and DE can be written as 


T m = T 0 and T d = T 0 ( 1 + uj d )exp 
So, we now obtain 


(2 + oj d ) \ da 


—3 I u> d + A , 

Ja 0 V (1 + u d )) a _ 


(43) 


A _ 3c 2 ( c 1 

Qm ~ 2 [ x ^r d n d + 


(44) 


and 


^ d ~~^r{ x W d +Ud+1 


^ h ~Hriw d +UJd] - 


(45) 


(46) 


Differentiating (45) and (46), we have 

Q d = lHcy/rr d \-^= { a(i + q-£=) + ^ i (1 - n d ) ( 1 + -y/n~ d ) - 3\n d \ U (a c 


-V^d 


y/rr d 3c 


V^d 


+ +11 1 + q) ~ 2 ( Wi ~ ')“ { (1 “ fij) l 1 + ■ 3Xitd 


(47) 


and 


Qh = c v^d |(1 ~ U d ) ^1 + ~ 3AD rf | — — \/Drf ) - 2(1 + ui d £l d )( 1 

+ qHRE ) 


(48) 
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FIG. 2. The plot shows the evolution of -pf- against A for three Planck data sets. 


Also, 


S E = t C 


( c 


4 Hvn~ d 


- 1 


(49) 


Then, from Eq. (25), we have 


dS% 3c 2 


dt 


1 


= — A-=- —+ 1 - 


27b V 


3c 2 


2?o(l + LOd) 


exp 


' ao 


+ A 


(2 + Ud) \ da 


(\ + UJd) J aj V yA^d 


A- 


+ + 1 


+ 


4 Hy/Tld Wn~d 
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- 1 + -AdH 


(yAh) ^ d ( ; 


yAh, 


+ ^d + 1 


-yAh, 


x { M1+9 ArJ + A^ 


| (1 ~ ^ d ) — 3Allrf| | + (A— j= + Wrf + 1) 


(4l 

l yAh, 


(l + ?)-2 




-1 - 


■\A2d 
2 


(1 — ^d) (1 + - \/^d) — 3A fid 


9 

- l AhH 


yAh 


x (1 + OJd^d) Cy/Qd | (1 — ^d) + -\AAj — 3AHd| — — y/VLdj — 2(1 + Ud^d) 


x (1 + qHRE ) 


(50) 


d.S$ 


Here again due to complicated mathematical expression for —rf-- we have shown its evolution 
against the parameter A graphically for three Planck data sets Q providing values for c and Qd- It 
is evident that GSLT holds for A > 0.8. 


IV. SUMMARY AND CONCLUSIONS 

The present work deals with non-equilibrium thermodynamics in the frame work of particle creation 
mechanism. At first, it has been shown that two fluids with different particle creation rates can be 
considered as a single fluid. However, if the dissipative pressure for the two fluids are identical in 
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magnitude but different in sign then the resulting single fluid has no particle creation rate and as 
a result it corresponds to a close system. Then considering the two dark species as DM (in the 
form of dust) and DE (as perfect fluid with a constant equation of state in one subsection and as 
holographic DE in another subsection), firstly we have shown that DM has a higher temperature at 
the initial stages and then it gradually decreases and attains a equilibrium era and finally the DE 
temperature dominates. Although the choice of the dissipative pressures (as in Eq. (29)) is purely 
phenomenological but from the mathematical point of view, it is easily solvable for both the choices 
of DE. Also physically we can argue as follows: As py oc H 1 2 3 4 5 6 7 8 9 10 11 , so dissipative pressure is of the order 


of H 2 and it is consistent with the choice of Barrow j2l|. Subsequently, in both the cases, analytic 
expressions for the total entropy variations have been determined. Due to complicated expressions we 
have graphically examined the validity of GSLT. In case of constant equation of state, GSLT holds for 
almost all values of A (except for very small A for ojy < — 1). For HDE model, we have chosen three 
Planck data sets for observed values of c and and it is observed that GSLT holds for all the three 
Planck data sets for A > 0.8. 
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